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^ 1. Introduction 
O 

^ , We study the properties of a particular one-parameter family of circle maps called off- 
(D , 

["t I I center reflections defined in §^. This map, in its 2-dimensional version, is first introduced in 

^ ' an open problem by S. T. Yau, 0, problem 21], who suggests a cross study of the dynamics 

, and geometry. In this article, we attempt to explore the possible link between the dynamics 

Q ' of this family of circle maps and their caustics. Although our study has not much contents 

I in differential geometry as Yau expected, it reveals some interesting phenomena. For example, 
we observe and partially prove that in a certain generic range of the parameter, the caustics 

I have exactly 4 cusp points for odd iterations; whereas for even iterations, each caustic is a curve 

o 

\G ' state about the dynamics and the geometry of the map; nevertheless, we still put it forward in 

o ■ 

0^ 



tangential to the circle at exactly four points. This may not be the best result that one could 



the hope that our study may invite better understanding to the subject. The off-center reflection 



I also bares several interesting analytic forms. It is a Blaschke product restricted to the circle. 

^ I Moreover, it has an infinite series expression highlighting that it is a perturbation of rotation on 

> 

X 



the circle. Since the work of Arnold, [A2], a standard type of perturbations has attracted much 



5-H ' interests in mathematics and physics communities, |BBJ, Di, 0]. This standard type is exactly 



a reduction of the series of the off-center reflection. This adds more flavor to our study. 

This family of off-center reflections plays an interesting role in the space of circle maps. 
With the parameter r going from to 1, it carries the initial antipodal map to the terminal 
doubling map, which provides a particularly nice way of deforming a simple dynamics to a 
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chaotic one (see the asymptotic orbit diagram in § |4.2| ). We expect that some notion of stabihty 
about the aforementioned cusp points and tangent points on the caustics might emerge from 
further study of this family of circle maps. We also hope to explore more systematically the 
dynamics of general circle maps through the method of symplectic and contact geometry in a 
forthcoming work. 

We will begin in §^ with the definition and some analytic properties of the map. Then 
the study is divided into three parts. In the first part our attention is given to the caustics 
of the map and its iterations. Our results concerning the caustics of odd iterations are more 
conclusive. Following a symplectic and contact geometry interpretation developed by Arnold 
I Al|, we discover the "generating function" for the corresponding Lagrangian embedding. This in 



turns provides explicit formulae for the orthotomics of the caustics. Moreover, these orthotomics 
are always smooth simple closed curves (and hence the caustic of an off-center reflection has at 
least 4 cusp points following the classical 4 vertex theorem), and they are convex in a certain 
range of the parameter. The method fails for even iterations. Nevertheless, explicit computations 
still provide reasonable support for our prediction. 

In the second part (^, the main observations are presented together with graphical 
illustrations. There are also theoretical support for them. For example, we have partial result 
that the caustics is stable with r < 1/3. The tedious computations are separated into §^]^ for 
detailed reading. 

In the third part (§P), we study the phenomenon of mode-locking behavior for this 
particular family of circle maps and the width of resonance zone is estimated. This is an 
attempt to understand the iterations of the map. This family extends a class of examples, which 
exhibits the same behavior, studied by Arnold and others. The mode-locking phenomenon of 
the off-center reflections and its "complex conjugates" are totally different. Moreover, r = 1/3 
is the first value that this behavior undergoes a structural change. This probably is not simply a 
coincidence with the bifurcation values of cusps. There is room for investigation in this aspect. 
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2. Off-Center Reflection 



An off-center reflection is a map §^ defined as follows: Pick a point, say (r, 0) in 



will be reflected at z with S as the curve of reflection and the reflected ray will hit another 
point Rr{z) on S-*^. The map z ^ Rr{z) : S-*^ ^ is what we call an off-center reflection. The 
action of the map is shown in the figure where a point on is represented by 4> mod 2it and 
a denotes the incident angle. 



The iteration of this map is a little uncommon at first sight because it is different from 
the usual successive reflection in a curved mirror. However, maps similar to it have been a 
center of discussion in circle dynamics. Let us first establish two analytic expressions for the 
map Rr :'S>^ ^ where < r < 1. By them, we may see Rr as a real function in terms of an 
infinite series as well as a complex function restricted to the circle. They are 



the interior of the unit disk . For any z G dD^ = S^, emit a ray from (r, 0) to z. The ray 




Rj.{(j)) = (j) + TT -2a 



mod 27r 



Q = a((/)) :^= Arg (cos (j) — r + \ sin (f)) — <j), 



where Arg is the principal argument taking values in (— vr, tt]. 



Lemma 2.1. The angle of incident a has the following Fourier sine series, 



oo 



.k 

— sin 



a((/)) = Arg (cos (f) — r + \ sin (p) — (p 




k=l 
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Thus, 



oo i. 



(1) Rr{(t)) = (l) + -K -2^ — sm(k(l)) mod 27r. 

fc=i 

Proof. It is clear that both a and da/dr are odd functions in (p. So, they both have Fourier 
sine series expansions. Let a = ^ ak{r) sm.{k(j)) . Then, the coefficients of the series for da/dr 
are given by 

(9ofc 2 sin0sin(A;(/)) 



dr TT Jq 1 — 2r cos </> + 

I r f cos{k - l)(f) cos{k + 



d(j) 

cnsf h -\- \ 
TT Jq \1 — 2r COS (j) + r^ 1 — 2r cos (p + r^ J 
1 / 7rr^~^ Trr^~^^ \ 



TT \1 — r^ 1 — r"^ ^ 
Hence, those coefficients for the series of a must be r^/k. □ 

In this paper, we often omit the mod 2tt when it is clear in context. Formula (|l|) of Rr 
without the modulo 2tt is exactly the lifting of Rr to a function from M to M taking to vr. 

By playing with the argument of a complex number, we get another expression for the 
map Rr- This shows that the off-center reflection is a special case of the so-called Blaschke 
product. This map is not extendible to the hyperbolic disk. 

Lemma 2.2. For \z\ < 1, we have a complex function 

2 1 - rz 



(2) -z 

z — r 

whose restriction to the boundary of the disk, \z\ = 1, is the off-center reflection Rj.. 

Therefore, this function Rr{4>) is harmonic when (r, 0) are treated as polar coordinates 
as it is the argument of the analytic function —z{l — z)^. We do not know whether this coupling 
between the parameters has additional physical or geometrical implications. 

By changing a sign of the off-center reflection, we have another map Rr defined by 

oo u 

— sm{k(l)) mod 2tt. 



k 

k=l 
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Geometrically, it is the mirror image of Rr{(j)) reflected by the diameter joining (p to (p + n. This 
map can be extended to the unit disk, namely, 

z — r 



1 — rz 

Therefore, it defines a map in PSL(2,M), the isometry group of the hyperbolic disk. The 
dynamics of Rr and this "conjugated" map Rr are completely different (see [p^ and §|5.2|). 



For sufficiently small r, Rr behaves very similarly to Rq, i.e. the antipodal map. In 
fact, when r < 1/3, it is in the same component of Rq in the group of orientation preserving 
diffeomorphisms of S^. However, i?i/3 is only a homeomorphism on and Rr a degree 1 map 
when 1/3 < r < 1. These can be easily concluded from the derivatives of Rr, which will also be 
useful later, 

1 — 4r cos (h + 3r^ 



R'ri^l^) 
Rri<P) 

Ri'\<P) 



1 — 2r cos (/> + ' 

2r(l — r^) sin (p 
(1 — 2r cos <j) + r'^y ' 

2r(l - r^) [(1 + r^) cos (j) - 2r(l + sin^ cj))] 



(1 — 2r cos (j) + r 



2\3 



In § [4.2| , we will give more information about the fixed point and other special points 
of Rr- More dynamical properties such as periodic cycles and whether they are attracting are 
discussed in |[Au[ . 

3. Caustics 

3.1. Caustic of Off-center Reflection. In this section, we will discuss the cusp phenomenon 
of the caustic of Rr- Classical examples of caustics are the locus of focal points with respect to a 
point on a surface and the focal curve of a convex plane curve. Corresponding to these caustics, 
there are the famous Geometric Theorem (Conjecture) of Jacobi and Four-vertex Theorem. 
There are many interesting at-least-four results, see ^ |A|, The caustic of off-center 



reflection provides another one. The conjugate locus of a point on a flat flying disc is, at 
degenerate situation, the caustic of the off-center reflection. It should be remarked that in the 



series of papers |BGG1, BGl, BG2, GK], the authors analyzed the singularities of the caustics 
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produced by a point light source when it is reflected in a codimension 1 "mirror" in and 
M^. Their emphasis though is on the "source genericity": whether the caustics could be made 
generic by moving the source. See also [BGG2]. 

For a circle map / : — > the family of lines joining cj) to /(</>) is 

F{(f),x,y) = (sin /((/)) — sin(/))(x — cost/)) — (cos /((/)) — cos </))(?/ — sincj)) 
= (sin /((/)) — sin (j))x — (cos /(</>) — cos (j))y — sin (/(</)) — (j)) ■ 

The caustic of the map / is defined to be the envelope of these lines. Thus, it is given by the 
dF 

equations -;—{(/), x,y) = = F{(j),x,y), that is, 



sin f {(p) — sin (j) — cos /((/>)+ cos \ fx 

/'(0) cos /((/)) -cos (/) /'((/)) sin /((/))- sin 0y \y 



sin(/(0)-</-) 
(/'(</.) -l)cos(/((/))-0) 



Solving for x,y, we obtain a parameterization of the caustic 

f'{(j)) cos (^ + cos /(0) 

X[(P) = 

(3) 



1 + /'('/') 
/'((/)) sin + sin f{(j)) 



1 + 

The tangent direction, which is degenerated at cusp points, of the caustic is given by 

/"(0)(cos - cos /(0)) - /'(</.)(! + /'((A)) (sin + sin /(</.)) 



(4) 



x'(c/>) 



(1 + /W)2 

/"(0)(sin - sin ficj))) + /' ((/<)(1 + nm^os cj) + cos /(</<)) 



(1 + /W)2 

The caustic (^ of the off-center reflection may run to infinity since 1 + R'r{4>) may be 
equal to zero. In fact, it is so if and only if r > 1/2. It should be more appropriate to define 
the caustic as the envelop of the geodesic normal field on the sphere. After the stereographic 
projection, it does not matter whether the caustic is defined on the plane or the sphere as the 
local properties remain unchanged (Darboux Theorem of symplectic structure). As we will see, 
the local property of the caustic of can be understood by direct computation. 



Theorem 3.1. For allO < r < 1, there are exactly 4 cusp points on the caustic of Rr- Two of 
them correspond to the Rr-orbit, {0, vr}, of period 2. 
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Proof. The derivatives of x and y can be expressed completely in terms of r and (p, namely, 

6r^(— cos (j) + r cos(2(/>))(r — cos (/)) siiKp 



(—1 — 2 + 3 r cos 0)^ 
6r^ (— 1 + 2r cos </)) (r — cos 0) sin^ i 



( — 1 — 2 + 3 r cos 0)^ 

The common solutions for x'{(j)) = = y'{(j)) are (p = 0,11 and two values of (j) with cos^ = r. 
Clearly, and vr are zeros of x' of first order and of y' of second order, thus, these are semicubical 
cusps. If cos(/> = r, after further differentiation and evaluation at the point, one has 
x» = -12r3; x"'(6)^'^^''"^''^- 



y"{<P) 



Qr\2r^ - 1) 



-6r3(10r2 - 3) 
1 — 



Thus, x"y"' — x"'y" = / 0. Therefore, there are also semicubical cusps at those values of (/) 
with cos 4> = r. □ 

Here are two pictures of the caustics of Rr, one for r < 1/2 and the other r > 1/2. Since the 
second one runs to infinity, it is drawn with "compressed" scale where a circle of radius > 1 
represents the point of infinity and the caustics has a self-intersection there. 



-1.5 



/ o.t 




\ o.t 


0.5 1 




caustic of Rr with r = 0.375 



caustic of Rr with r ~ 1/ 



3.2. Symplectic and contact geometry interpretation. The explicit computation in the 
previous section has its advantages and shortcomings. On the one hand, it gives a very exact 
count of the number of cusp points; on the other hand, it is too specific and also very complicated 
to apply especially when iterations are considered. In this section, we present the symplectic 
method which may work in general, though the result obtained is not as specific as before. The 



terminology of [Al] is followed. 
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Denote the coordinates of the unit cotangent bundle ST 



by (p 

X ; Py 5 y ) where 



{x,y) G and + = 1- This bundle is a contact 3-manifold with the contact 1-form 
Pxdx + Pydy, and the cotangent manifold T*(M^) is symplectic with the symplectic 2- form 
d{jpx dx+pydy). 

With our notation of (j) and a in §|2[ the vector from (p to Rr{4>) is never zero and the 
unit vector in that direction {px,Py) is given by 

Px = cos{(p + vr — a), Py = sm{(j) + vr — a). 

Then 

fpx\ I cos(0 + 7r-a) \ 
Py _ sin(</) + vr — a) 

X cos (j) + S cos{4> + Ti — a) 

\y ) \ sin (j) + S sin(</) + vr — a) y 

defines a map L : x — > T*{M?). This map may be thought of as a flow (in parameter S) 
of unit speed in the direction of the reflection lines, starting with the round circle 5" = 0. This 



is a case of what Arnold called "Legendrian collapsing" |A1|. 

Let p : r*(R^) — > be the canonical projection. The Jacobian of the map p o L is 
J(p o L) — det ' + vr — a) — sin (p — S{1 — a') sm{(j) + vr — a) 



sin((/) + IT — a) cos (j) + S{1 — a') cos{(j) + vr — a) 
cos{(j) + vr — a) cos cp + sm{(j) + vr — a) sin p + S{1 — a') 

— cos a + S{1 — a') 



So the equation for the critical curve on §^ x M} is 



S 



cos a 
1-a' 



Proposition 3.2. The critical curve, when mapped to the {x,y)-plane, agrees with the caustic 
of Rr. 



Proof. Notice that since RriP) = </> + vr — 2a, we have 

R'^ = l-2a' or 1 + R'^. = 2(1 - a'). 
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Therefore, 

, cos a , , , 

X = cos (p H cos(0 + vr — a j 

1 — a' 

(1 + i?^) cos + 2 cos a cos(0 + vr — a) 

~ i + i?; 

_ cos + cos((/) + vr — 2a) 

which agrees with one of the equations (Q). We may get the expression for y similarly. □ 

Furthermore, a straightforward calculation shows that 

■px dx + pydy = sin ad(j) + dS. 

Thus, the image of L is a Lagrange cylinder in T*(R^). Notice further that a = a(0) is an odd 
function of (j) and therefore 

/ sinadcf) = 0. 

This implies that p{L) is an exact Lagrange cylinder. Take a function S{(j)) along the circle by 

r<t' 

S{4>) = — sin a d(j). 



Jo 

It defines a section in the Lagrange cylinder p{L). It is easy to see that —£'((/>) is increasing for 
< < vr and decreasing when vr < < 27r. Therefore, the curve C given by 

X = cos (p + S{(j)) cos(0 + TT — a) 

y = sin 4> + S{4>) sin((/) + vr — a) 

is quite likely to be a convex plane curve. Let us show that this is the case when r < 1/2. 

First we note that C has a continuous normal field (cos(0 + vr — a), sin(0 + vr — a)). It 
is easy to compute 

x'^ + y'"^ = 4sin2 a + (cos a - S{1 - a')f . 

Therefore, x'^ + = is possible only when </> = vr. But the number of zeros of x'^ + should 
be even (geometrically, because C is co-oriented). Thus x'^ + y'^ > all the time and C is a 
smooth simple closed curve. We also point out that the curvature of C is 

1-a' 
^Jx'^ + y'2 
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which is non-negative when r < 1/2 and therefore C is convex. Now the family of reflection lines 
of Rr is identical to the family of normal lines of this convex curve C. Therefore, the caustic of 



Rr has at least 4 cusp points [Al]. 

The function S{(p) should be thought of as the generating function of the circle map 
Rr{(j))- The curve C is related to the orthotomic of such a reflection. We will study such 
generating functions for general circle maps in a forthcoming work. 

3.3. Iterations of Reflections. For an integer n, we denote the iteration of Rr by ii" = 
Rr o i?"^^ : ^ S^. The equations (^) and (^) in §3T, with / = R^, give a parametrization of 
the caustic of i?" and its tangent. 

The cusps on caustics of ii" are more intriging and complicated than that of Rr- There 
are fundamental differences between the caustics when n is odd or even. To see this difference, 
we may consider the trivial example that r = 0. For any odd n, Rq is the antipodal map and 
its caustic is a point (a degenerated curve with cusp). However, for even n, we have the identity 
map, so the caustic is defined by the family of tangents and it is the circle itself (a smooth 
curve). It is expected that this cusp versus "smooth" situation remains for r close to 0. 

Theorem 3.3. For sufficiently small r > 0, the caustic of i?^™"*"^ has at least 4 cusp points. 



Proof. To some extent, the symplectic method in § p.2| may be adopted for i?^'""'"^. It can be 
proved by induction on m that 

R'r+\<P) = (l) + 7r-2am{<P) 

for some odd function «„,((/)). For example, 

ai{(l)) = a((/)) + a{(j) + vr - 2a{4>)) + a{(f) - 2a{4>) - 2a{(f) + vr - 2a{(f>))). 

Therefore, we still have an exact Lagrangian cylinder and a sectional curve defined by 

f't' 

S{4') = — sind,; 
Jo 



'■m ' 



For small r, it is a convex curve and hence there are at least 4 cusp points on the caustic. □ 
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Remark. The argument fails for even iterations of Rr, because the analog of a is not an odd 
function. Furthermore, it is not clear about how small the range of r should be. Yet, from 
experimental observation, there are at least four cusps for any r > and there are exactly four 



for < r < 1/3. For more, please see the discussion after proposition 4.3. 



4. Experiments and Observations 

To get more accurate information about the caustics of iterations, we have to rely on 
lengthy calculations. Our investigation is indeed partly theoretical and partly experimental. 
We will first describe some interesting properties with illustrations. The technical details of 



justification are left to the interested reader in §4.2 



4.1. Observations. We first put forward a conjecture about the exact picture of the caustic 
when i?" is still a diffeomorphism. Then we look at the bifurcation process of the structure of 
the caustics when r varies. Finally, we compare the caustics for different n. We will soon see 



that the 2-cycles of Rr play a special role (propositions 4^ and 4^). The 2-cycles are {0, vr} 



1 _ ^1 + 8r2 

and { ib(/)c } where (pc G (0, vr) and cos (pc = . We will often refer to this notation. 



Conjecture 4.1. For < r < 1/3, the caustic of Rr"^~^^ is a CP° curve with exactly four cusp 
singularities, with two of them occurring at (f) = 0, vr. On the other hand, the caustic of R?/"" is 
a differentiahle curve; CP" everywhere except at exactly the four 2-periodic points of Rr, where 
the caustic is tangent to the unit circle. 



The conjecture is demonstrated by the following pictures, which are produced by pro- 
gramming in Mathematica. The purple (thin) curve is the unit circle and the blue (thick) curve 
is the caustic with the light source at the red dot, (r, 0). Besides these pictures, we also have 
partial results that support our conjecture. 



12 



T. K. K. AU AND X. S. LIN 




caustic for n — 2, r — 1/3. 



caustic for n — 3, r — 1/3. 



caustic for n — 4, r — 1/3. 



Proposition 4.2. Any caustic curve of R'^'^ is tangential to the unit circle at (cos sin at 
any point cp satisfying = cp. In particular, this includes the points 0, ir, and ±(pc- 

Moreover, if r < 1/3, these are the only four points that the caustic meets the circle. 

Proof. By substitution of such (j) into equations @ with / = i?", we have x{(f>) = cos cj) and 



y{(f>) = sincj). Applying lemma 4^ to such (p's, the assertion about the tangential property of 
the caustics of even iterations of Rr follows easily. Secondly, since + = 1 , we have 

(RrYicI)) cos (/) - cos R'^ic/)) = 
sine/)- sin = 0. 



This leads to R^{(j)) = 4>. If r < 1/3, by lemma 4.8, n must be even and (f) is one of the four 
2-periodic points. □ 

We have already seen from symplectic topology that caustics of odd iterations, R"^."^^^, always 
have at least four cusps for sufficiently small r > 0. Moreover, for all Rr, two of the cusps occur 
at (/> = 0, vr. Now, we may extend this result to odd iterations of Rr with isolated exceptional 
values of r. It turns out these exceptional values only occur in r > 1/3. 

Proposition 4.3. For < r < 1/3 and for generic 1/3 < r < 1, the caustic of R^"^'^^ always 
has cusps at 0,TT. 

Proof. Note that for both (pa = 0, vr, one has Rr{4>a) = 4>a mod 2tt, thus we may apply 
lemma 4.5 to check whether there are cusps. Furthermore, Rr{(j)a) = (pa and R['{(pa) = 0, we 
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can simplify the chain rules (lemma |4.6| ) and obtain {R'^^'^^)" {(j)a) = and 

{R^+^)^'HM = K{Rr{M){K)'{<Pa)B^\<t^a) + 

+ K{ct>af [{R':)'{(l>a)Rf\Rr{<Pa)) + KiRr{<Pa}f (R'^Y'HM 

Let us temporarily define, for positive integers n, 

A^ = -{R^y + {R^y^ + 2{R'^)^^l 

So it suffices to check that A2m+i{4'a) 7^ 0. We will proceed by induction on m. First, by direct 
computation, 

MO) = (Y37j2 > 0' M^) = > 0. 

Then, it can be shown that 

- R'Ma)R'r{<Pa + ^) + R'AM^ R'r{<Pa + T^f] [RrY {<Pa) 



1 - 9r' 

I — f2 



2\ 3 



AM+A2{<Pa)-{R'^y{4>a). 



In particular. 



^2m+l(0) = ^ ^2m-l(0) +^2(0) . ^ ; > 1 ^ 

'1 — r^y \1 — rj\l + r 

l-9r2\\ .^ , , , /I - 3r\"'"Vl + 3r 



^2m+l(^) = 2" ^2m-l(vr) + A2(7r) 

'1 — r^/ Vl — \l + r 



where 



[1 — r'^y 

It can be easily computed that ^2(0) > for < r < 1/3 and A2{t^) > for all r > 0. Thus, 

1_ 9^2x3 



A2m+l{(t>a) > yY~i^ J "^2m-l(0a) > 0. 

For r > 1/3, vl2m+i may have zeros. We can only conclude from above that ^2m+i is a 
rational function in r with denominator being a power of (1 — r^). Therefore, it has only isolated 
zeros and cusps at 0, vr occur for r > 1/3 generically. □ 
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We observe that at cusp points, @ gives a set of "homogeneous" equations which has zero 
determinant. Thus, except at a couple of 0's, it is sufficient to solve only one equation, say, 
x'((j)) = 0. It is likely that this equation has exactly four solutions for < r < 1/3. However, 
it is still hard to solve explicitly especially for high iterations of Rr- Here are two pictures of 
the functions {Rf)"/{{Rf.y{l + (Rf)')) and (sin</. + sin i?3((/)))/(cos (/> - cosi?3(0)) for r = 0.1 
and r = 0.33. Exactly four cusp solutions are demonstrated in each of them. We do not have 
a proof of this graphical fact. Perhaps it may be proved by detailed curve sketching argument 
and comparison of i?^™""^ and i?^™"^^ using the known properties of R^ given in lemma |4.8| . 





It is probably worthwhile to compare the situation here with the so-called Jacobi conjecture 



promoted by Arnold [Al, A5|. As mentioned in the beginning of §3J-, the caustics of Rr agrees 
with the locus of conjugate points of (r, 0) on a flat flying disk. Although the loci of higher order 
conjugate points are not the same as the caustics of odd iterations of Rr, their common contact 
geometric nature indicates that two problems of whether there are exactly 4 cusps on the loci 
of higher order conjugate points and whether there are exactly 4 cusps of the caustics of odd 
iteration of Rr might be related. In both cases, since we do not have exact nice formulae for the 
loci of higher order conjugate points and the caustics of odd iterations of Rr, it would be very 
difficult to have an exact count of cusps. Of course. Conjecture [4.1| deals with a very special 
situation. From graphical evidence, it is tempting to think that a proof should not be out of 
reach with brute force calculation. Nevertheless, after spending much effort on this temptation, 
we think some more conceptual understanding of the caustics of iterations of Rr is needed in 
order to get an exact count of cusps. 
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Furthermore, in the proof above, we see that for n > 3, except A^^n), the quantity 
^n(O) = = An^n) always at r = 1/3. There is a possible structural change on the caustic 
of i?" occurring at r = 1/3 for n > 3. We observe from experiment that, for the caustic of 
odd iterations, once r > 1/3, bifurcation of cusp may occur. Interestingly, from the computed 
pictures (only that of iij? is shown), bifurcation only occur at the cusp corresponding to (p = 
but not others. Would the different properties between ^2(0) and ^2(71") be part of the reasons? 




975 -0.95 



n — 3, r — 0.36. its enlargement, — 0. 

On the other hand, bifurcation into cusps also occurs for even iterations of Rr at r = 1 /3. 
We will discuss this by beginning with some Taylor expansions. Since Rr and its iterations are 
27r-periodic odd functions, they have particular nice expansions at (pa = 0, vr. This enables us to 
see the local properties of the caustics more clearly. 

Let / be any even iteration of Rr, then f{(pa) = 4>a- We write ^ = (j) — 4>a and g{'d) = 
f{(t)) — (pa and suppose it has an expansion 



2k+l 



gi.^) = ^a2fc+ii9 

k=0 

One can inductively work out the coefficients of the expansions of l+f'{(j)) = 1 +(/'(??), cos /(</>) 
lb cosg{'d), etc. If P^, Qk denote polynomials with Pfc(0, ... , 0) = = Qk{0, ■ ■ ■ j 0), one has 



1 + ai (1 + air + 



+ ai ^ [l + aiY 
These expansions are helpful to understand the caustics of i?^™ at (pa = 0, vr. It would be 
convenient to look at the pictures before we go on. 
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In the above pictures, cusps are born near cj) = and vr. From the enlargement, the caustic 
bifurcates into 2m cusps when r increases across 1/3, where Rr changes from a diffeomorphism 

] . By this, we wih see that r = 1/3 
is the value that the caustic of i?^™ changes at (p = 0, it. In fact, the caustics of Ry^, 
following Taylor expansions. At (j)a = 0, 

x((/)) - x(0) = ^1?^ + 0(i?^) 
?/(</')- 2/(0) = 18# + 0(i?5); 

and at (j)a = IT, 

243 

x(0)-x(7r) = — 1?^ + 0(^6) 
lb 

y(</.)-y(^) = ^# + 0(^5). 

This shows that the caustic of R^ undergoes a swallowtail bifurcation at 0, vr when r = 1/3. We 
may further work out the expansion of R'fJl as the m- iteration of Ry^- Using oi = and 03 7^ 
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for R^, we have 

for some function U with U (0) 7^ 0. The bifurcation of the caustics of i?^™" at (/) = 0, vr should 
be of the type (t?^'""*"^, ■i?'^'") when r passes 1/3. 

In the above, we consider the behavior of the caustics with parameter r and n fixed. 
What happens to the caustics if r is fixed and n is allowed to vary? There is an interesting 
phenomenon for r < 1/3. Although there is fundamental difference when n is odd and even, 
this difference disappears as n goes to infinity. Here are the pictures showing how the caustics 
of and i?^™ change-they tend to the same quadrilateral. 




caustics of -R°/3 caustics of 



The green (dashed) vertical line from (pc to —(pc in the picture of even or odd caustics is deter- 
mined by r but not the number of iterations. When n is even, it is where the caustic is tangent 
to the circle. When n is odd, every caustic is tangent to this vertical line because it is the line 
jommg (j)c and = The point of tan gency occurs exactly at (j)c by definition. 

Proposition 4.4. For < r < 1/3, as m ^ 00, both the caustics of R'^."^^^ and i?^™ approach 
the same quadrilateral defined by the four points 0, vr, and i(j)c, which are the only 2-periodic 
points of Rr ■ 

Proof. Let us first see that the caustics of i?^™""*"^ at the four points tend to the circle as m — > 00. 
These four (p's are the solution to Rl{(t)) = (/). Thus, we have = i2;(</.)™+i • 

RURr{(l>)r. 
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At (/> = 0, TT, the coordinates of the caustic are 

It is clear that (i?2™+i)'(0) = R'M""^^ K{-^T = (l^)™^^ • (w)™- '^^^S' ^(0) ^ 
m ^ cxD. The situation at vr is similar. 

At the point </>c with i?"((/>c) = — </>c, the coordinates of the caustic at this point are 

X(<^e) = COS C^., = ^ ^ . 

Since i?r. is odd and R'^{4'c) = (pc it follows that {R^)'{(l)c) = R'j.{4>c)^. Moreover, by 4rcosi;^c = 



1 — vT+~8r^, one may show that 

(-l + (i?-y(0e)) 

1 + {R^yi^c) 

Again, these two cusps approach to the unit circle. 



1 as n — > CO. 



Secondly, from the lemmas 4^ and 4^, and vr are the attracting fixed points of R^ while 
ib(/)c are repelling. Moreover, the attracting basins for and vr are (— 0ci ^Pc) and ((/)c, 2vr — </)c) 
respectively. Thus, for any given neighborhood of 0, for sufficiently large m, for any neighboring 
4'1j4'2 G {—4'c,4'c), R^^{(t^i) and i?^™((/»2) lie in that neighborhood of 0. Hence, the intersection 
of the lines from to R'^'^{(j)j) lies in a neighborhood of the quadrilateral. The proof for the 
cases at vr and of odd iterations are similar. □ 

4.2. Technical Results. In this section, a couple of technical results will be given. They are 
mostly done by direct computations and the methods may not be insightful. However, they may 



be the necessary evil for they will be used to justify our observations in § |4.1| . The first one deals 
with the existence of cusp at certain special "symmetric" positions. 

Lemma 4.5. Let f denote any iteration of R^. On the caustic of f , the conditions for the 
occurrence of a semicuhical cusp at (j) are 

. /'(</,) = and/"(0)/O i//(0) = 0; 

. = and -f'i^) + + 2/(3) (0) / if f{<P) = + vr; 
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Proof. This is proved by computing the derivatives of (|3|) and (Q) , then evaluate at the particular 
values (pQ or (pa- Then the result follows by verifying x' = = y' and x"y"' — x"'y" 7^ 0. □ 

Remarks. Although this lemma is stated for an iteration of Rr, it is actually true for any circle 
map. Likewise, many results in this section hold in a more general setting but we would like to 
focus on iterations of Rr- Furthermore, at a point (/) with /{(/)) = ~4': we always have x'{4>) = 0. 
The conditions are 

f'mi + /'(0)) cos ct> + r(0) sin = 

f'{(p){2 + cos(2(/))) + 6f'{(t>f cos^ + (1 + 2 cos{2(t>)) f {(t>f - 2f^^\ct>) sin^ / 0. 

Analogously, if f{(t)) = n — (p, we have y'{4') = and conditions 

/'(<A)(1 + /'(</.)) sin <A -/"(</.) cos (^ = 
/'(0)(2 - cos(20)) + 6f'{^f sin^ + (1 - 2 cos{2^)) f {^f - 2f^\<t)) cos^ / 0. 

The next one may be an exercise for calculus students. 

Lemma 4.6. The chain rules for i?" with n = p + q are given by 

(R^Yicp) = {RpyiR'^m ■ {Riy{(P) = r'm ■ KiRrm K{K-\<t^)) 

(Kn^) = iR?)"iRK<p)) ■ iRiri<pf + m'iR'rm ■ mn^) 
(Ri^f^H^) = (RpfHR^m ■ iR'rn<pf + {R^r)'{Rim ■ m^''\<t>) + 
+?,m"{Rim-{Ri)'{<p)-{Ri)"{<p). 

In determining the cusps on the caustic, some orbits in the iteration play a special role. We 
thus establish the following to handle that. 

Lemma 4.7. For < r < 1/3, i?^™'^^ has no fixed point and / + tt. 



Proof. Firstly, one can obtain algebraically the four fixed points of R^. The attracting ones are 
0, vr, while ^(pc are repelling. Then by simple calculus, we get the corresponding attracting 
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basins and conclude that i?^"*((/>) converges to or vr mod 27r monotonically. By the series 
expression (^) of iJ^ one can deduce the following estimate 

\Rr{(^)-(l)\>7r-2\logil-r)\. 

Then, according to the convergence of -R^"^ ((/)), the lemma follows. □ 



/ 1 _ ^1 + 8r2 

Lemma 4.8. For < r < 1/3, let (pc = arccos ^ |, the solution sets to the 



4r 



following equations are given as, 



• Rr{4>) = (j) has solutions if and only if n is even, which are 0,7r,ib(/)c. 

• Rr{(f)) = —(j) has solutions ziz(j)c if n is odd; and 0, vr if n is even. 

• Rri(f>) = (p + TT has solutions if and only if n is odd, two of them are and tt. 

• R^{(j)) = TT — (p has solutions 0,7r when n is odd and two solutions when n is even. 

Proof. First of all, the preceding lemma already give us part of the conclusions. Furthermore, 
R^ fixes the intervals [— vr, —(l)c\, [— ^oO], [0, (/)c], and [(/>c,7r] in a way that 

Rl{(t)) < (p on (-vr, -0c) or (0, (j)c), 
Rl{<j)) > (j) on (-(/)c,0) or ((/fcvr). 

Thus, induction process like R^^{4)) < R'^'^"'^ (ip) < (p on the corresponding intervals will give 
the only fixed points of -R,^™. The other claims are done similarly. □ 



The above lemmas are illustrated by the following figures. 




Plots of Rl, Rl and R^^, r = 1/3. 



1 2 3 4 5 

Plots of Rl, R^ and R^, r = 1/3. 
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As an attempt to understand more about the iterations i?", we computed the asymptotic 
orbits of the two critical points it arccos (^^-^^f—^ for 1/3 < r < 1. The asymptotic orbit of (p 
is the set { R^{(p) : A'^i < n < A''2 } for large Ni,N2. The plot of asymptotic orbits against the 
parameter r is called a bifurcation diagram. It shows the attracting periodic cycles or chaotic 
behavior of a map according to the variation of the parameter. These diagrams have been 



further analyzed in [ Au 
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5. Mode-locking 

5.1. Background. The study of circle maps is closely related to the study of differential equa- 
tions on torus (i.e., equations with double periodic coefficients). For any such differential equa- 
tion, one may consider the Poincare return map of the flow, which defines a map on a meridian 
circle of the torus. It turns out that the stability of the equation is reflected by this circle map. 

Far back in 1959, in his Ph.D. project, Arnold investigated the circle map 

4> ^ (j) + a + e cos 6 



and obtained information on its resonance zone in the (a, e)-plane, ||A2| , §12]. This gives rise to 
the famous picture of so-called Arnold tongues. Subsequently, there are numerous studies, by 
physicists and mathematicians, [BBJ, Di, FKP, IBB, |P|, on the perturbation of a rotation 



(/!> 1-^ (/) + — e sin </), ee[0, 1). 
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The focus is on the phenomenon called mode-locking and the Devil's staircase. Arnold later 



gave a proof in [A4] of his observation for circle maps of the form 

4> (j) + 0, + e (trigonometric polynomial) 



as well as analytic reduction of many circle maps, [A3, Ch. 3, §12]. The algebraic nature of the 
method is also apparent in the problem of particular differential equations. Arnold predicts that 
a general theorem exists for these equations and general circle maps. 

In this section, we will provide further evidence towards Arnold's prediction by showing 
similar behavior in the off-center reflection. It should be remarked that the off-center reflection 
is not of the form studied by Arnold. Thus, it may be another small step towards the general 
theory. 

We consider a two-parameter model of circle maps which arises from the off-center re- 
flection map, namely, with parameters r £ [0, 1) and J7 G (— 7r,7r], 



OO ^ 



RrM^I') =(t> + n-2Y^ — sm{k(l)). 

k=l 

Here we use r instead of e to be consistent with previous sections. Note that unlike the models 
discussed above, r cannot be factored out. This map can be thought of as an imperfect off-center 
reflection on the circle where the reflected angle has a constant deviation from the incident angle. 
It is the original off-center reflection when 17 = vr. We may not get such a deviation by varying 
the metric of the circle; it is better understood in terms of symplectic geometry. 

For (pQ G S"*^, there is the rotation number uj{Rr n, 'Po) = li™ — ~ ! where the 

' n^oo n 

right hand side is performed on a lifting of Rr^n- It is independent of 0o if Rr,n is diffeomorphic. 
In such case, one simply denotes Lu{Rr^n). If Rr,n is only a degree 1 map, one has a rotation 
interval instead. These notions are indeed defined for any circle map. Historically, attention 
has been centred around perturbations of rotations, i— > + J7 + u{(j)). It is natural to ask for 
the relation between and uj. The physicists usually refer to Q as internal frequency and to 
as resonance frequency. When uj = uj{^}) is a locally constant function, the situation is called 
mode-locking. 



OFF-CENTER REFLECTIONS 23 

Herman has extensively studied the mode-locking property and obtained interesting re- 



sults, I HI, H2]. These results are applicable to Rr^n because it satisfies the property Aq defined 



by Herman. 

Theorem 5.1. For all ujq G 27rQ and < r < 1/3, there is an interval I = 2r of loq such that 
for every O G Ir, the diffeomorphism Rr^n has rotation number ujq. 

The interval I is called resonance interval and its size depends on r (and of course 
iOo). Its variance in terms of r defines a picture which looks like a tongue. We will discuss it 
later. Furthermore, from Herman's study, the off-center reflection model also demonstrates the 
well-known Devil's staircase. 

Theorem 5.2. For any < r < 1/3, the function Q i— > cj(i?r,rj) is nondecreasing, locally 
constant at any rational number, and has a Cantor set of discontinuity. 

We have mentioned that if we alter a sign and form the "conjugate" family 

k=l 

the dynamics is completely different. Actually, Rrn can be extended to e^'^'^ on the 

1 — rz 

hyperbolic disk, which defines a hyperbolic element in PSL(2, M). The mode-locking phenomenon 
does not occur, i.e., u! [Rr,n) = 'ipir/q only if 17 = p/q. 



5.2. Width of Resonance Zone. In |A4|, Arnold discusses the mode-locking situation of a 



rotation slightly perturbed by a trigonometric polynomial, g{x), 

f : X ^ X + + eg{x). 

The resonance zone is the set {(J7,e) : Q GXg}. Arnold developed a formal calculation to 
estimate the width of the interval in terms of e, which gives rises to a picture of the resonance 
zone. This formal calculation is related to the homological equation of analytical reduction. 



|A3|. If the rotation number is rational, the width of the resonance interval 2^ is bounded by a 
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power of £. The graphical plot of the resonance zone in the eO-plane form the so-called Arnold's 
tongue. 

By a method similar to Arnold's, one may also estimate the width of for the off-center 

reflections Rr,n, < r < 1/3. We will show the different behaviors of R = Rr^-n and R = i?r,7r at 
the same time. 

For simplicity of computation, let us first consider the resonance zone containing tt. 
Writing = tt + a, the second iterates of the maps are 

oo ^ oo ^ y k \ 

R^[x) =x + 27r + 2a-2^^sin(fcx) -2^^sinfc f x + 7r + a - 2^ ^ sin(fex) j 

k=i k=i ^ ^ 

oo I, CO I, / b \ 

R (x) = x + 27r + 2a + 2 sin(fcx) + 22j-^sinfe(x + 7r + a + 22_]-^ sin(fcx) j 

k=l k=l ^ ^ 

The equations of resonance are R? {x) = x + 2tt and Pl' {x) = x + 2'k. Let v = a^2Y^'^ sin(A;a:) , 
we have 

oo 
k=l 

where v = vir + V2r'^ + v^r^ H . Note that the solutions of v's for R and R do not only differ 



— sin(A;2:) — sin k(x + v) 

k k 



by a sign. One can see this by the subtle combinations of the signs of the infinite series in their 

— 2 

second iterates. Inductively, one may show that for i?, we have '^k = -j^ sin(A;x), while those for 
R are 

Vl = — 2sin(a;), 
V2 = sin(2x), 

g 

vs = 2sin(x) — - sin(3a;). 



The above result leads to a r-series for a and its maximum and minimum provide bounds for 

the resonance zone, namely, 

{2sin(2x)r^ + [2sin(x) — | sin(3x)] + ■ ■ ■ , for the map R, 
for the map R. 

This calculation agrees with our previous remark that mode-locking (near lo = ir) does not occur 
for R. Furthermore, we have 
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Theorem 5.3. The width ofir is bounded by Cr^ for J7 = tt and Cr for general 

The computation for the resonance zone at a general O = 2p7r/q is more comphcated. 
The equation to formally expand is -^^ a+2p7r/q(^) = ^ + ^Ptt. The coefficients of 

a = air + 027-^ + a^r^ + • • • 
provide the estimates of Ir- It turns out that the first term ai does not vanish, indeed, 

q-l 



qai = 2 y sm I x H 1 



j=0 

This may not be a sharp estimate, yet we can only conclude that the width of Xj. is of order r 
in general. 
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